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l> : ABSTRACT 

ly-^ ■ We present the results of a simulation to investigate the prospects of measur- 

^ ■ ing mass, age, radius, metallicity and luminosity data for brown dwarfs in fully 

Psj ■ eclipsing binary systems around dwarf spectral types from late K to early M that 

could be identified by ultra-wide-field transit surveys such as SuperWASP. These 
O ■ surveys will monitor approximately a million K and M dwarfs with |6| > 20° 

(where blending of sources is not a significant problem) at a level sufficient to 
detect transits of low luminosity companions. We look at the current observa- 
tional evidence for such systems, and suggest that ~1% of late K and early-mid M 
c/3 ■ dwarfs could have a very close (~0.02AU) BD companion. With this assumption, 

and using SuperWASP as an example, our simulation predicts that ~400 brown 
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dwarfs in fully eclipsing binary systems could be discovered. All of these eclips- 
ing binaries could yield accurate brown dwarf mass and radius measurements 
from radial velocity and follow-up light curve measurements. By inferring the 
brown dwarf effective temperature distribution, assuming a uniform age spread 
and an a=0.5 companion brown dwarf mass function, the simulation estimates 
that brown dwarf brightness could also be measurable (at the 10% level) for ~60 
of these binary systems from near infrared follow-up light curves of the secondary 
eclipse. We consider irradiation of these brown dwarfs by their primary stars, 
and conclude that it would be below the 10% level for ~70% of them. This means 
that in these cases the measured brown dwarf brightnesses should essentially be 
the same as those of free-floating counterparts. The predicted age distribution of 
the primaries is dominated by young systems, and ~20 binaries could be younger 
than IGyr. Irradiation will be below the 10% level for ~80% of these. We suggest 
that many of these young binary systems will be members of "kinematic moving 
groups" , allowing their ages to be accurately constrained. 
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Subject headings: stars: low-mass, brown dwarfs - stars: fundamental parameters 
- binaries: eclipsing - techniques: photometric 

1. Introduction 

The mass of a brown dwarf (BD) is arguably its most basic characteristic. The essence of 
the classification of an object as a BD depends on the objects sustaining deuterium burning 
(for solar metaUicity >0.013Mo) but not sustaining hydrogen burning (for solar metaUicity 
<O.O75M0). However, since BDs contract, cool and fade with time, their observable proper- 
ties depend strongly on both age and mass. Furthermore, because BD colours are strongly 
influenced by atmospheric dust formation and molecular absorbers (like H2O), their appear- 
ance will depend not only on effective temperature (Tg//) and surface gravity (g), but also 
on composition. Accurate BD mass, age, radius, metaUicity and luminosity data are thus 
vital to properly test our understanding of these extreme objects. 

In view of the particular importance of such data for BDs, it is unfortunate that it has 
thus far not been possible to completely characterise any BDs in this way. There are several 
reasons for this. In order to measure BD masses, one must find BDs in binary systems, and 
measure the orbits. One also needs to be able to measure the BD itself to determine its 
luminosity, and the age of the binary must be established, so as to infer the age of the BD. 
Proper motion measurements have shown that ~10% of solar type stars have BD companions 
at wide (>1000AU) separations (Gizis et al. 2001). However, the periods of such binaries 
are much too long to be useful in measuring masses. A small number of possible close BD 
companions to solar type stars have been identified (Halbwachs et al. 2000; Nidever et al. 
2002; Blundell et al. 2004). But such systems arc very rare, and could not be directly 
imaged anyway because of the glare of the primary. For the lowest mass stars, adaptive 
optics imaging has shown that close (~1AU) BD companions to late M and early L dwarfs 
are significantly more common (Close et al. 2003), with orbits that can be measured in 
a few years, and binary components that can be individually studied. For example, Bouy 
et al. (2004) astrometrically measured 36% of the 10 yr orbit of 2MASS 0746+2000AB, 
constraining the system mass at the 10% level. However, the ages of such binaries will in 
general be very difficult to constrain without recourse to the theoretical models that one 
wishes to test. For instance, Bouy et al. compared 2MASS 0746+2000 A and B to Dusty 
model predictions (Chabrier et al. 2000) in a K, J-K colour-magnitude diagram (CMD) and 
estimated a system age of 150-500Myrs. However, when compared to late M and L dwarfs 
with accurate parallaxes, 2MASS 0746-1-2000 A and B do not he above the CMD sequence 
as sources of this age should. A potentially more useful type of system is represented by the 
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Epsilon Indi multiple, which consists of a double BD binary (visual separation of 2.65AU) 
at a distance of 1500AU from its K5 primary (McCaughrean et al. 2004). The BD binary 
orbit should be measurable (P~15yrs), and the age of the system can be constrained to 
some degree (0.8-2Gyrs) by measurements of the K dwarf primary (Lachaume et al. 1999). 
However, such multiple systems are not common place, and age constraints established from 
measurements of a single primary star will not generally be very precise. 

The concept presented here is to exploit the new breed of ultra- wide-field transit surveys 
(designed primarily to find signatures of planetary transits) to identify BDs in fully eclipsing 
close binary systems around cool stars. BD masses could be measured from the dopplcr 
wobble of the primary (using radial velocity techniques). One could measure BD radii 
from the depth of the primary eclipse (when the BD crosses the face of the primary). The 
brightness of the BD may be measured from the secondary eclipse (when the BD goes behind 
the primary) provided the BD is bright enough compared to the primary. Many of these 
systems will be young (because BDs are brighter when young), making membership of young 
kinematic moving groups more likely. Membership of such moving groups would accurately 
constrain age and metallicity. 

In this paper, we discuss the prospects of identifying BDs in binaries capable of providing 
the full compliment of BD mass, age, radius, metallicity and luminosity data. In Section 2 we 
discuss what is currently known about the frequency and separation of close BD companions, 
and how this effects the likelihood of eclipse. In Section 3 we use the SuperWASP survey as 
an example of the new ultra-widc-ficld transit surveys, estimate sensitivities and the effects of 
source crowding in images with very large pixels, and estimate the distances out to which BD 
transits should be detectable. Section 4 describes the additional factors limiting the number 
of suitable binary systems; the accuracy with which secondary eclipses can be measured 
(Section 4.1), and the limits imposed by the doppler wobble technique (Section 4.2). Section 
5 gives some background on kinematic moving groups, and details how one could accurately 
constrain the ages of young systems. Section 6 describes our simulation to predict the size 
and properties of the eclipsing BD population. Section 7 then discusses some identification 
and foUowup procedures, and Section 8 presents our conclusions. 

2. Close BD companions 

2.1. Frequency and separation 

In Table 1 we show and compare what is known about the frequency of companion 
stars and BDs as a function of separation. For solar type stars, Halbwachs et al. (2003) 
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and Duquennoy h Mayor (1991) have measured the frequency of stellar companions as a 
function of period (which we convert into separation), using 10 years of radial velocity data 
combined with visual binaries and common proper motion systems. It is clear from this work 
that stellar companions are common over a large range of separation. However, it has been 
known for some time that there is a dearth of low mass ratio (q<0.1) companions in close 
orbits around solar type stars, known as the "BD desert" . Other radial velocity surveys (eg. 
Marcy & Benitz 1989; Fischer & Marcy 1992; Nidever et al. 2002) have now studied large 
numbers of solar type stars, and it is currently estimated that < 1/500 have a BD companion 
within 3AU (Vogt et al. 2002). This corresponds to a BD companion fraction of <0.07% per 
decade in separation (clearly much less than the ~4% for stellar companions. The desert 
does not extend to wide separations, with sub-stellar companions imaged at 15-20AU (Els 
et al. 2001; Lui et al. 2002), and BDs are as common in very wide binaries (~1000AU or 
more) as stellar companions around solar type stars. 

In terms of primary type, there is evidence (from the OGLE transit survey) that the 
desert extends out to the mid-K dwarfs, since the majority of the small low luminosity com- 
panions discovered by Udalski et al. (2002a, b, 2003) turned out to be low-mass stars and 
not BDs (Dreizler et al. 2002). However, for later dwarfs, the target lists for radial veloc- 
ity surveys are not large enough to accurately constrain the close BD companion fraction. 
Despite these limitations some interesting candidate BD companions have been discovered. 
Reid & Mahoney (2000) found a possible BD companion (0.07-0.095Mq) at a separation of 
~0. 005-0. 007AU from an ~M5 dwarf in a radial velocity survey of ~50 Hyades M dwarfs. 
Nidever et al. (2002) found 2 possible BD companions (at separations of ~0.2 and ~2AU) 
around early M dwarfs from a radial velocity survey including ~100 M dwarfs. Conversely, 
Marchal (2003) found no close BD companions in a comparable radial velocity survey of 
~100 MO-5 dwarfs. 

At somewhat larger separations, Oppenheimer et al. (2001) found 1 BD in a corona- 
graphic survey of the stellar population within 8 parsecs of the Sun. This volume limited 
sample (ie. mostly M dwarfs) consisted of 163 stars in 111 star systems, and BDs with 
mass>O.O4M0 and age <5Gyrs could have been detected with separations of 40-120AU 
around 80% of the sample. BDs are not found in significantly wider orbits around early M 
dwarfs, perhaps because impulsive perturbations by close stellar encounters will cause such 
weakly bound systems to dissolve (see Burgasser et al. 2003). It is not, then, surprising 
that this lack of wide BD companions is also seen around the latest stars and BDs (with the 
possible exception of one very young candidate BD binary which would have a separation of 
~240AU; Luhman 2004). However, closer companions to late M, L and T dwarfs are quite 
common (see Table 2), and while accurate masses are generally not available, many of these 
will be sub-stellar (eg. Freed et al. 2003; an L7.5 BD ~3AU from an M8 star). 
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Of interest here are the late K-mid M dwarf primaries, since we will subsequently show 
that sufficient numbers of these dwarfs will be monitored by ultra-wide-field surveys to find 
eclipsing BD companions (see Section 6). It is clear from Tabic 1 that BD companions 
to these stars could occur at the ~2% level in orbits close enough to make transits likely 
(<0.02AU; see Section 2.2). However, an accurate determination of the masses of the BD 
candidates discovered so far would require knowledge of orbital inclination. Some may turn 
out to be low-mass stars at low inchnation (eg. Halbwachs et al. 2000). If correct though, 
it would show that although the BD companion fraction in close orbits around such stars is 
lower than that of stellar companions, there would be no BD desert for these primaries. 

In addition to these systematic searches, there have been two serendipitous discoveries 
of possible close BD companions. Schuh et al. (2003) reported the discovery of an echpsing 
BD companion to a late K dwarf with a separation of 0.02AU. This system was discovered 
while photometrically monitoring a white dwarf. Also Santos et al. (2002) discovered a 
radial velocity BD companion to an M2 dwarf with a separation of 0.017AU. This system 
was found as a result of the weak gravitational wobble it induced in the visual K2 companion 
of the M dwarf. 

Having considered all known candidates, it can be seen that the spectral type range of 
the primaries in systems that could potentially harbour significant populations of very close 
BD companions ranges from late K to mid M. The orbital periods of these BD companions 
in close orbits around late K-mid M dwarfs is 1-1.5 days. At these separations, a full transit 
lasts ~1.5hrs. 



2.2. Eclipsing fraction 

The majority of the known candidate BD companions around K and M dwarfs (three out 
of five) are in close orbits with separations <0.02AU (see previous section). Close systems 
such as these will be the most likely to eclipse, because eclipsing probability is oc 1/a (where a 
is separation). We therefore chose to consider the likelihood of an eclipse for a BD companion 
at a separation of 0.02AU. This separation corresponds to ~10R* for an early M dwarf, and 
for such separations, a favourable ~1 in 15 should be eclipsing. The requirement for a full 
eclipse is that 

sin(90-i)< ^^~^^^ , (1) 

where i is the inclination angle, Rp and Rbd are the primary and BD radii respectively, and 
a the component separation. The fraction of sources that should be eclipsing is (90 — i)/90. 
Using theoretical stellar radii from the 500Myr NextGen models of Baraffe et al. (1998), and 
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assuming O.IR© for the BDs (typical for 500Myr age), we get the echpse fractions shown 
in Figure 1. As the primary radius decreases, the chances of a full eclipse fall off, and for 
comparable radii full eclipses are very unlikely. It can be seen that for late K-early M dwarfs, 
the eclipsing fraction is ~4-7%, and so even with a ~1% binary fraction at 0.02AU, one could 
still expect to find fully eclipsing BD companions around ~1 in 2000 dwarfs of this type. 

2.3. Irradiation 

Since we are interested in measuring companion BD properties that would be essen- 
tially the same as those of free-floating counterparts, we desire irradiation effects to be low. 
However, for the very close BD companions discussed in Section 2.2, this may not always 
be the case, since the cool BD companion could be heated by the nearby brighter hotter 
primary, as happens with giant planets in close orbits around solar type stars (hot Jupiters). 
Hot Jupiters are strongly irradiated, and have higher Te//S and swollen radii as a result 
(eg. Burrows Sudarsky & Hubbard 2003; Chabrier et al. 2004). Irradiation should only be 
an issue when the incident flux from the primary is significant compared to the emergent 
flux from the companion. If the incident flux is 10% of the emergent flux for instance, then 
the BD luminosity would increase by no more than 10% (and the radius and Teff by much 
less than this). This level represents a useful level of accuracy with which to measure BD 
brightness. Figure 2 shows how the amount of irradiation varies for BD companions (spectral 
types M6-T4) at 0.02AU separation from their primary stars. Each point in the plotting 
area represents a binary where the primary and secondary spectral types can be read off 
the y and x axis respectively. Using NextGen model data for the primaries, BD Te//s ap- 
propriate to their spectral type (Reid et al. 1999; Knapp 2004) and assumed BD radii of 
O.IR0, we derived a locus of binary systems in which the BD will be irradiated by 10% of its 
emergent flux. This locus is plotted in the flgure as a solid line, where binaries above the line 
are not significantly irradiated. For example, an M8 BD will not be significantly irradiated 
at 0.02AU from a mid K dwarf (or later) primary, and an L4 BD will not be significantly 
irradiated at 0.02AU from an MO-2 (or later) primary. 

2.4. Age and mass distributions 

The mass and age distribution of the companion BDs is important because the older 
and lower mass BDs will be more difficult to measure from the secondary eclipse. The age 
distribution should be the same as the local disk. An analysis of chromosphcric activity of G 
and M dwarfs in the solar vicinity (Soderblom, Duncan & Johnson 1991; Henry, Soderblom & 
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Donahue 1996) is consistent with a uniform star formation rate over the age of the disk, and 
we thus assume that BDs will have a uniform age distribution between and lOGyrs. The 
mass distribution of companion BDs is less constrained. The BD mass function slope for the 
field and young clusters/associations ranges from a ~0.5-l (where dn/dm oc m~"). However, 
the masses of BDs in close binaries may have a very different distribution, since dynamical 
factors during formation (eg. Reipurth & Clarke 2001) and possible radial migration (due 
to early disk interaction; Armitage and Bonnell 2002) could be important. For very low- 
mass binary systems for instance, the q distribution peaks at 0.7 - 1 (Burgasser et al. 2003; 
Pinfield et al. 2003). Such a bias amongst binaries with an early M dwarf primary would 
produce a strong preference for stellar companions over BDs, and the BD companion mass 
function itself would prefer high mass BDs. However, it is not currently possible to estimate 
a for close BD companions from the limited observational evidence, and it essentially remains 
a free parameter. Low values of a would result in proportionally more high mass BDs (for 
some value of the BD companion fraction), which would be brighter and easier to detect. 
Higher values would mean proportionally more low-mass BDs, many of which would be too 
faint to detect. 



3. Ultra- wide- field transit surveys 

The discovery of the first transiting extra-solar planet (Charbonneau et al. 2000; Henry 
et al. 2000) has led to a proliferation of ground based photometric searches aiming to detect 
planets around bright solar type stars by searching for transit signatures in photometric light 
curves (see Horne 2003). These searches fall into two categories; the deeper, but relatively 
small area searches such as the Optical Gravitational Lensing Experiment (which uses a 1.3m 
aperture telescope), and the ultra-wide- field surveys such as Super WASP (Christian et al. 
2004), TrES (Alonso et al. 2004), Vulcan South (Caldwell et al. 2003) and others, that use 
small aperture (~10 cm) wide-field (5-10°) CCD based systems with spatial resolution of 
~10" pixel"^. While accurate light curve extraction is more challenging for ultra- wide-field 
data, the ultra-wide-field surveys promise to discover transits around brighter nearer stars, 
for which close transiting companions can be more accurately measured. In this paper we 
focus on the abihties of these ultra-wide-field surveys, and use the UK's SuperWASP project 
as a specific example. 



-8- 



3.1. SuperWASP 

SuperWASP (Super Wide Angle Search for Planets) is a project to construct and op- 
erate facilities to carry out ultra-wide-field unfiltered photometric monitoring surveys (see 
PoUacco et al, www.superwasp.org). The first SuperWASP facility (SuperWASP-I) is lo- 
cated on Roque de Los Muchachos, La Palma, and a second facility has also been funded 
for the southern hemisphere. Super WASP-I consists of five camera units mounted on a large 
computer controlled fork mounting, housed in a small enclosure. Each camera unit consists 
of a 200mm fl.8 Canon lens with an aperture of 11.1cm (made of fiuorite and UD glass com- 
ponents) and a 2048x2048 pixel e2v42 thinned CCD (BV coated) cooled to -60C by a 3-stage 
Peltier. Each camera has an 8°x8° field of view, resulting in large 14.3 arcsecond CCD pixels. 
The combined sky coverage of the five cameras is thus ~320 sq. degs. SuperWASP-I current 
observing strategy is to repeatedly measure 30s integration images of the same region of sky 
throughout each night for a period of ~2 months, before moving onto a new region. With 
low overheads (~4s readout). Super WASP-I takes ~100 sets of images per hour. First hght 
was in November 2003, and at time of writing Super WASP-I is routinely taking data. A 
data reduction pipeline is being implemented, and light curves of all stars measured will be 
placed on the public "Ledas database" at The University of Leicester, where software search 
algorithms will be available for data mining. 



3.2. SuperWASP sensitivities 

SuperWASP is an unfiltered system, and as such the transmission profile (or pass band) 
is determined by the atmosphere, the optical throughput of the lenses and CCD entrance 
window, and the CCD quantum efficiency. To estimate the shape of the SuperWASP trans- 
mission profile we have used available site and technical details. We used La Palma atmo- 
spheric extinctions from the "INT Wide Field Survey" web-site (based on observations by 
Derek Jones in August 1993). For the Canon lenses, the optical transmission of UD glass and 
fiuorite (~95% per pair of air-glass interfaces) does not vary significantly from ~0.23-7 mi- 
crons (minimising chromospheric aberration). We thus assumed a wavelength independent 
optical transmission of 75% to allow for the lens groups (in each Canon lens) and the CCD 
entrance window. The CCD quantum efficiency was taken from the E2V CCD selection 
guide, interpolating between the -20°C and -100°C curves. Figure 3 shows the atmospheric, 
lens and CCD transmission profiles (dotted lines) against wavelength, and the resulting to- 
tal SuperWASP transmission (solid line). The transmission profile peaks in the V and R 
bands (at ~60%). Blue of V, both atmospheric extinction and the CCD quantum efficiency 
reduce throughput. Red of R the transmission is predominantly reduced by CCD quantum 
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efEciency. We estimated a system zero point using a spectrum of Vega (from Hayes, 1985 
and Mountain et al., 1985), which we take to represent a zero magnitude source. Accounting 
for the collecting area (96.7 sq. cm for each lens) and the Super WASP transmission, we in- 
tegrated the detected Vega photons across the SuperWASP band. The predicted zero point 
was found to be 21 magnitudes (2.5x10^ detected counts/s from a zero magnitude source). 
As a guide, this means that a signal-to-noise of ~50 would result from a 30s integration on 
a 16*'* magnitude source. 

3.3. Unfiltered colour term 

We used the spectra of Pickles (1985; available on CDS) and Leggett et al. (2000; 
from www.jach.hawaii.edu/~skl), combined with a Johnson V-band (Bessell 1990) and the 
SuperWASP pass band, to derived synthetic V-band and SuperWASP-band flux ratios for 
early G-mid M dwarfs. These flux ratios were converted into SuperWASP colours (V-m^^y, 
where m^^ is a magnitude on the natural SuperWASP system) using our Vega spectrum to 
define zero colour. Figure 4 shows the resulting colour terms as a function of spectral type. 
It can be seen that for early G dwarfs the colour term is small, but for later spectral types it 
becomes significant because most of the optical flux from these dwarfs is not in the V-band, 
but is in the red part of the SuperWASP band. 

3.4. Unfiltered sky brightness 

To estimate the La Palma sky brightness as seen by SuperWASP, we used the full 
range of optical sky brightness measurements from Benn & Ellison (1998), for bright, grey 
and dark sky conditions. For each sky we converted the magnitudes into fluxes using our 
Vega spectrum (see Section 3.1). The resulting flux levels were then joined together to 
provide an approximate sky spectrum. These sky spectra were then multiplied through by 
the SuperWASP band pass, and converted into SuperWASP magnitudes using our system 
zero point. The SuperWASP sky brightness estimates were found to be 18.3, 20.3 and 21.1 
magnitudes per square arcsecond for bright, grey and dark skies respectively. If one accounts 
for the large SuperWASP pixels and a spatial resolution element containing 12 pixels (a 2 
pixel radius aperture), one would expect sky brightnesses of 9.8, 11.8 and 12.6 magnitudes 
per resolution element, for bright, grey and dark skies respectively. 
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3.5. Source crowding with large pixels 



The large SuperWASP pixels (and consequent large aperture sizes used for photometry) 
mean that image crowding will be much more of an issue than it is for higher spatial resolution 
imaging, particularly in the Galactic plane. This is important for two reasons. Firstly, if 
a transit source shares its aperture with another star (ie. if it is blended) , the SNR of the 
transit detection will be decreased. However, more importantly, if the PSF of a neighbouring 
source spills into the aperture of a star, this fraction of the measured flux could vary over 
time if the PSF changes slightly, or the chosen aperture does not remain perfectly centred 
on the star - SuperWASP tracking (~0.01arcsecond/s accuracy) will cause stars to drift 
between pixels over several hours. Since we are interested in identifying eclipse dips that 
are ~4% deep, we would not want such flux pollution much above the 1% level. We have 
therefore used synthetic gaussian point-spread-functions to determined the distance at which 
a source (with FWHM=1.5 pixels) will pollute it's neighbours 2 pixel aperture by 1%, and 
estabhshed how this distance {di%) depends on the relative brightness of the two sources 
(Amag). We found that 



for Amag <2.0. For fainter neighbours, di% drops more rapidly, reaching zero for Amag —5 
(ie. a factor of 100). So if a source has a neighbour with d < di%, we consider it significantly 
blended, and not measured with sufficient accuracy for our purposes by SuperWASP. We 
applied this criteria to a series of 1 sq. deg catalogues (centred on different Galactic coordi- 
nates) from the SuperCOSMOS Sky Survey database^. Since we are interested in late K and 
early M dwarfs, we only considered SuperCOSMOS sources with Bj-R>1.8, and estimated 
the fraction of blends as a function of Galactic coordinate and source brightness. We found 
that for V-16, -30% of sources with \b\ > 20°, and ~75% of sources with |6| < 20° will 
be blended in SuperWASP images. For V~18, we found 50% and 85% blended fractions for 
|6| > and |6| < 20° respectively. 



di% = FWHM X (1.95 - 0.12Ama^), 



(2) 



3.6. Transit detection magnitude limits 



The signal-to-noise level of a transit detection is 





^http: / /www- wfau.roe.ac.uk/sss/ 
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where Rbd and are the radii of the companion BD and the primary star, jS'obs and noiscobs 
are the signal detected from the somT:c and the associated noise during an observation, P is 
orbital period, tobs is the time taken for an observation (~30s), and a is the separation of 
the binary. The first term represents the relative decrease in brightness during transit, the 
second term is the signal-to-noise associated with each observation, and the last term is the 
square root of the number of observations made during a transit. The noise model we used 
was 

noiseobs = ^J (-^obsiobs) + (-^sky^obs) + (^^obs) + (NR^) + (Syst^), (4) 

where ^sky is the signal detected from the sky in an N=12 (2 pixel radius) spatial resolution 
element (of. the lenses produce a ~1.5 pixel point-spread- function) during an observation, 
D is the dark current (~0.01e~/s), R is the read-noise (~12e~), and Syst represents any 
systematic errors. In general, systematic uncertainties associated with relative photometry 
of bright sources can be reduced to the 1-2 m-mag level (Kane et al. 2004; Bakos et al. 
2004). For the fainter sources, of interest here, systematics will be dominated by field 
crowding effects which sets the upper level of the systematic uncertainties at a self imposed 
10 m-mags (see Section 3.4). 

We define our transit detection limits requiring a SNR>5. We would only expect 1 in 
~2 million stars to register a 5-a transit due to random noise alone, and since SuperWASP 
will survey rather less than this number of late K and M dwarfs, false positives from random 
noise will not be a problem. We combined equations 2 and 3 with the SuperWASP zero point, 
and used the binary characteristics and stellar model data described in Section 2.2 to derive 
SuperWASP magnitude hmits for 5-a transit detection. These limits, for bright, grey and 
dark skies are shown in Figure 5 (solid lines), as a function of primary My- Primary spectral 
types are also indicated at their appropriate position (taken from Pickles 1985). V-band 
magnitude limits were then derived using the colour term from Section 3.2 (dashed lines). 
For comparison, the G star lO-cx bright sky detection limit of Horne (2003) is over-plotted 
(open circle), and joined to its 5-cr equivalent (filled circle). It is clear from the figure that 
as one moves to later type stars, transits can be detected out to fainter magnitudes. This is 
because later stars are smaller, and eclipse dips arc correspondingly deeper. Also, the effect 
of the broad unfiltered pass band is clear, since SuperWASP will do better for red sources 
than a comparable V-band filtered system. 

3.7. Distance limits for transit detections 

The V-band transit detection limits from Figure 5 were converted directly into distance 
limits using the appropriate values of My. These distance limits are shown in Figure 6 (solid 
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lines) for bright, grey and dark skies. It can be seen that transits of late K-early M dwarfs 
should be detectable out to ~200-500pc. For the later mid M dwarfs, SuperWASP is limited 
to ~100pc, due to their fainter intrinsic brightness. The figure also contains some additional 
limits, which will be described in the next section. 

4. Limits imposed by follow up requirements 

4.1. Measuring BD brightness from secondary eclipses 

Having identified a good transit signature of a low luminosity fully eclipsing companion 
(and confirmed it with an intermediate sized optical telescope) , one can then attempt to 
measure the brightness of the companion from the secondary eclipse. Our ability to do this is 
governed by the accuracy with which we can measure follow up light curves and the relative 
brightness of the companion and primary. Compared to late K/early M dwarfs, BDs with 
spectral types <T2 are relatively brightest in the K-band (see Leggett et al. 2002). K is 
generally thus the best band in which to search for a secondary eclipse. With My=8-12 
and Mk=5-7 (see Reid & Cruz 2002), a typical late K/early M dwarf near the SuperWASP 
distance limits (500-200pc respectively) would have K^13.5. The accuracy of K-band light 
curves will depend on photon noise, stability of the instrumental response, and atmospheric 
conditions. Using 2-4m class telescopes, a signal-to-photon-noise ratio of ~1000 could be 
achievable for a K~13.5 source in ~5-10 minutes. Variations in the zero point and extinction 
can be corrected for using frequent standard measurements, flat-fielding should be good to 
at least 0.1% accuracy, and atmospheric scintillation noise for such observations will be ~1 
part in 10^ (equation 10 of Dravins et al. 1998). The main source of systematic uncertainty 
in near infrared photometry is the amount of atmospheric H2O absorption, which can be 
variable on very short timescales. Kidgcr & Martin-Luis (2003) measured high precision near 
infrared photometry of bright (K<10) northern hemisphere stars using a photometer, and 
were able to obtain K-band accuracies (inferred from multiple measurements) as good as 
0.4m-mags for some sources (see their Figure 8). Typical uncertainties would not always be 
this good, as atmospheric conditions change (cf. Kidger & Martin-Luis' median uncertainty 
was ~3m-mags), however an overall K-band measurement accuracy of Im-mag should be a 
reahstic level to aim for. Note that large scale NIR arrays (such as WFCAM on the UK 
IR Telescope) would allow comparison stars to be measured at the same time as targets, 
which could reduce atmospheric systematics. However, we will continue to use Kidger & 
Martin-Luis as our main example, since they have observationally established these levels of 
accuracy. 

The total accuracy with which we can measure the depth of a secondary eclipse dip in a 
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near infrared light curve (atot) will depend on both the accuracy of an individual observation 
{(Jobs) a-nd on the number of data points (N) measured during secondary eclipse (where 
o'tot = (JobsVN ; see Henry et al. 2000). For the binaries we are considering, we would expect 
an eclipse to last ~1.5 hrs, and we could thus expect to obtain 6-10 K-band measurements 
per eclipse (allowing for 5 minutes overhead per observation). By measuring 3-4 eclipses, 
and phasing the data appropriately, one could then measure the depth of the secondary 
eclipse dip with an accuracy of 0.2m-mags. Thus, for an eclipse dip to be measured with 
10% accuracy, it must be at least 2m-mags deep. The BD must therefore be no less than 
500 times (6.75 mags) fainter than the primary star. Estimating primary from Reid & 
Cruz (2002), and BD Mk (appropriate for M6-T4 spectral types) from Dahn et al. (2002) 
and Knapp et al. (2004), we have derived My brightness limits for primaries with Mk 6.75 
magnitudes brighter than BD companions of each spectral type. These primary My limits 
are shown in Figure 6, where each is labelled with its appropriate BD spectral type. As 
an example, it should be possible to measure the secondary eclipse of an L6 eclipsing BD 
companion of an My >9 primary. 



4.2. Measuring BD masses from radial velocities 

If a low luminosity eclipsing companion is detectable from its secondary eclipse, then 
radial velocity techniques can be employed to accurately measure its mass. For eclipsing 
systems, the inclination is known to better than 1°, and the relative mass of the BD (com- 
pared to that of the primary) can be determined by measuring a radial velocity curve for the 
primary. Primary masses will thus act as anchors for the BD masses, and uncertainties in 
the mass-luminosity relation for low-mass stars (cf. Henry & McCarthy 1993) will propagate 
through to the BD masses. However, it is expected that the low-mass star mass-luminosity 
relation should improve significantly in the next few years, as the orbits of many more M 
dwarf systems are measured (eg. by the RECONS; Jao et al. 2003). 

The radial velocity accuracy required to determine BD masses places its own limits on 
the distance of useful binary systems. To illustrate this, a O.O5M0 brown dwarf at 0.02AU 
will cause a radial velocity oscillation of ±20kms~^ for a O.5M0 M dwarf primary. Standard 
arc lamp echelle spectroscopy is limited by several factors; zero point drifts in the wavelength 
scale between exposures, the spectral features available for cross correlation techniques, the 
presence of telluric absorption, the signal-to-noise of the spectra themselves, and intrinsic 
properties of the star (surface spots and/or rotation). We will return to the issue of spots 
and rotation in Section 7.3. For the other factors however, we will use previous optical 
studies with 8m telescopes as a guide. Reid & Mahoney (2000) observed zero point drifts 
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of ±lkms~^ during a full night for the Keck HIRES echelle spectrograph. However, they 
managed to obtain ±300ms~^ accuracy by cross correlating echelle orders (containing strong 
telluric signatures), with white dwarf spectra (which act essentially as continuum sources 
onto which a clear telluric signature is imprinted). More recently, Konacki et al. (2003) 
used an iodine cell to more rapidly measure zero point drifts of ~200ms^^ over timescales of 
~30 minutes. These studies looked mostly at K and M dwarfs from 3850-6200A and 6350- 
8730A respectively, and found that such drifts generally dominate the overall uncertainty, 
with contributions of only ±60ms~^ coming from the signal-to-noise (~15-25) of the spectra 
themselves. 

It is thus clear that optical echelle spectroscopy should be capable of providing radial 
velocities with accuracies of ±0.5kms~^ (enough to constrain close BD companion masses at 
the ~2-3% level) using lower signal-to-noise spectra from rather fainter sources with V~19. 
We have converted this magnitude limit into a distance limit which is shown in Figure 6 
(dashed line). It can be seen that for sources with spectral types <M4 the echelle followup 
requirements do not impose any additional distance constraints beyond those of SuperWASP 
itself. The efficient follow up of later type primaries at the SuperWASP distance limit would 
require a near infrared echelle spectrograph (such as PHOENIX on Gemini). However, we 
do not expect to find many of these as fully eclipsing systems (see Sections 2.2 and 6). 



5. Determining accurate ages 

5.1. Open clusters 

Disk stars with the best constrained ages arc those in open clusters. The age of an open 
cluster can be well constrained using a variety of methods, such as fitting the upper main 
sequence turn off (Sarajedini et al.l999), or (for young clusters) measuring the magnitude of 
the lithium depletion edge (Stauffer, Schultz & Kirkpatrick 1998). An open cluster is beheved 
to be a coeval association (with only a relatively small age spread of perhaps a few Myrs). 
Therefore, if membership of a cluster is established, then the cluster age can be assumed 
for the star. However, the majority of disk stars are not open cluster members. Within the 
SuperWASP distance limits (~500pc) there are ~31 open clusters in the La Palma sky (see 
the LEDAS database^) . Although catalogued membership of these clusters is by no means 
complete, a reasonable estimate would be that between them they contain no more than a 
few thousand late K/early M dwarfs. This represents a tiny fraction of the number of non 



^http: / /ledas- www. star .le.ac.uk 
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cluster disk stars of this type in the Super WASP sky (~1 miUion). Open cluster membership 
will thus not be very useful for constraining the age of a significant fraction of the binary 
systems under consideration. 

5.2. Young moving groups 

Open clusters are thought to represent only small components of larger kinematically 
distinct coeval populations that are spatially dispersed within the disk. Such disk populations 
arc known as moving groups. Moving groups were originally identified by Eggen in the 1950s. 
However, it has since become clear that many of the older (>lGyr) so called moving groups 
he identified show a spread in both age and mctallicity (Nordstrom et al. 2004), and these 
kinematic signatures arc more likely the result of disk heating by stochastic spiral waves (eg. 
De Simone, Wu & Tremaine 2004). There is evidence however, that the younger moving 
groups are coeval. Castro, Porto de Mello & de Silva (1999) found identical Cu and Ba 
abundances for 7 G and K dwarfs in the Ursa Major moving group, and King et al. (2003) 
found clear evidence of a correlation between membership of the Ursa Major moving group 
and activity levels (from Call H and K emission strength). 

Such young moving groups are thought to originate in the same environments as open 
clusters. Very young clusters should virialise somewhat within their natal gas, and when the 
gas is cleared by OB star action and the cluster expands, the stars with high enough velocities 
will become unbound. These unbound stars then slowly expand forming the moving group, 
and possibly leaving behind an open cluster consisting of any remaining bound stars (see the 
N-body simulations of Kroupa Aarseth & Hurley 2001). Typical expansion velocities for the 
unbound population will be ~5-10kms~^ (5-lOpc/Myr), and the moving group space motion 
will thus remain quite well defined (superimposed on the more dispersed disk population) 
until it has had time to become significantly phase mixed by disk heating. Models suggest 
that phase mixing will increase a moving group velocity dispersion by ~20kms~^) over ~lGyr 
(Asian, Figueras & Torra 1999). Moving groups younger than ~lGyr should thus consist of 
young populations with characteristic space motions, and membership of such a group will 
accurately constrain the age of an eclipsing BD-M dwarf binary just as it would do for open 
cluster members. 

The most notable such groups in the solar neighbourhood are hsted in Table 2 along 
with their kinematic motions and age. Most have open cluster(s) associated with them, 
which range in size from the rich Pleiades cluster down to the almost evaporated Ursa-Major 
cluster. In general, each group has one (or two in the case of the Pleiades) distinct kinematic 
signature. However, note that the Gould belt is rather more complex than the others, since 
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it consists of an agglomeration of several very young kinematic groups superposed on a less 
distinct velocity ellipsoid (Moreno, Alfaro & Franco, 1999). Each of the distinct kinematic 
sub-groups is believed to have a common origin (eg. Asian et al., 1999), with the ellipsoid 
possibly resulting from an expanding ring of star formation (the local bubble). Interestingly, 
it appears that a large fraction of young local stars reside in these moving groups. Dehnen 
(1998) has shown (using Hipparcos astrometry) that the majority of colour selected young 
disk stars are moving group members, and Makarov (2003) found that a high fraction of 
X-ray active stars within 50pc belong to the Pleiades and Gould Belt moving groups. 



5.3. Identifying young potential moving group members 

In order to separate possible young moving group members from the large population 
of older disk stars, one needs to observationally identify young sources. Age constraints can 
be placed on young late K and early M dwarfs in a variety of ways. Chromospheric/coronal 
Ra emission from open cluster members (Hawley et al. 1999; Gizis, Reid & Hawley 2002) 
have shown that there is a well defined, age dependent colour beyond which activity be- 
comes ubiquitous. Despite considerable scatter, all stars redder than this colour are dMe 
(EWijQ, > l.oA), while the bluer stars are dM without emission. This so called "Hq, limit" 
colour increases with the age of the population, and can identify (in emission) K5-M0 dwarfs 
younger than ~30Myrs, MO-3 dwarfs younger than lOOMyrs, and M3-5 dwarfs younger than 
~600Myrs. Lithium abundances can also provide age constraints. Late K dwarfs deplete 
lithium over ~100Myrs, and early-mid M dwarfs deplete lithium over ~30-50Myrs (Preibisch 
et al. 2001; Zapatero-Osorio et al. 2002; Jeffries et al. 2003; Randich et al 2000). One may 
also constrain ages by fitting isochrones to colour- magnitude data. A ^O.SM© star will 
not have fully contracted until it is ~100Myrs old. Therefore, with an accurate parallax 
distance (possible even at 500pc using the VLT as an interferometer), an appropriate colour- 
magnitude diagram (see StaufFer et al. 2003), and photometric corrections to account for 
metallicity sensitive colours (eg. Kotoneva, Flynn & Jimenez 2002), the age of such a star 
could be approximately constrained. 

Clearly, identifying late K and early M dwarfs younger than ~100Myrs will not be 
difficult. However, older counterparts will be fully contracted, show no lithium and have 
H„ emission <lA. Such sources will require higher resolution analysis if their ages are to 
be constrained from their activity levels. For instance, one could gauge activity levels by 
measuring Mg II h & k and Ca II H and K (Doyle et al. 1994; Christian et al 2001; 
Mathioudakis & Doyle 1991), or by measuring at high resolution, and look for weaker 
line emission in the core of the absorption line (eg. Short & Doyle 1998). These methods 
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allow one to study activity levels down to those associated with old disk M dwarfs (see Doyle 
et al. 1994), and should therefore provide an effective way to separate young (<lGyr) late 
K and early M dwarfs from older disk stars. 

5.4. Kinematic membership of moving groups 

Establishing kinematic membership of moving groups will require space motions accu- 
rate to a few kms~^. Since accurate radial velocity curves will have been previously measured 

(to determine BD mass; see Section 4.2), these will provide very accurate centre of mass ra- 
dial velocities for the binaries. Although it should be possible to derive some proper motions 
(with the required accuracy) from existing archive measurements (eg. with SuperCOSMOS) 
sources with distances near 500pc would require some additional effort. Proper motion ac- 
curacies would need to be at the milli-arcsecond/year level (cf. 5kms^^ ~2mas/yr at 500pc), 
and it is now quite feasible to measure proper motions of this accuracy using adaptive optics 
technology with a baseline of ~l-2 years (cf. 3mas positional accuracy; Close et al.., 2003). 

The confidence with which one can assign kinematic membership of a particular moving 
group will depend on how unique the group's measured kinematic signature is compared to 
other young moving groups and any non moving group population of young stars (see Figures 
1 and 2 of Montes et al. 2001). The Hyades moving group for instance has a space motion 
well separated from the local standard of rest, but has a rather larger internal dispersion 
compared to some younger moving groups. As a result some Hyades moving group members 
will not be obvious kinematic members of the group. At the other extreme, the very young 
members of the Gould belt and the Pleiades moving group components will not show as 
much internal dispersion (within each component). However, due to their similar space 
motions, there is some kinematic overlap between some of these components. Furthermore, 
young stars that are not members of any moving group can act as a source of potential 
contamination, if they have the same space motion as the group. However, if the majority 
of young stars are moving group members (see Section 5.2), such contamination should not 
be at a high level, and will only decrease membership probabilities by a small amount. 

While it is clear that only some fraction of moving group members will be easily identi- 
fiable from their kinematics, it is by no means clear exactly what this fraction will be, since 
there are insufficient radial velocity measurements of young disk stars to properly charac- 
terise our expectations. However, it is worth noting that moving group characterisation will 
certainly improve in the near future, as large scale radial velocity surveys (such as RaVE; 
Steinmetz 2003) provide radial velocity and metallicity for complete magnitude limited sam- 
ples. 
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6. Simulating the number of useful eclipsing BD companions 

We will now describe a simulation we have made to predict the size and properties of 
the population of fully eclipsing BD companions that could be identified by the SuperWASP 
survey over the whole sky. This simulation predicts how many BDs will be identified, how 
many of these will be detectable from their secondary eclipse, how many are likely to be 
significantly irradiated by their primary, and how many will be young (and thus potentially 
moving group members). 

6.1. The simulation 

For our simulation, we consider low-mass stars with BD companions in close (~0.02AU) 
orbits. These binaries will have short orbital periods of 1-1. 5d. We thus expect ~40-60 
transits/binary during a ~2 month SuperWASP monitoring period, ~11-16 of which should 
occur while SuperWASP is taking data (allowing for typical 9h nights and ~25% time lost to 
bad weather). SuperWASP should experience ~9.5 dark nights per lunation, and the moon 
will be set for typically ~20% of non-dark nights. We therefore expect dark sky conditions 
~45% of the time, and would expect ~5-7 transits to occur during this time. It is therefore 
appropriate to use the dark sky distance limits for transit detection. We also assumed that 
for late K dwarfs (V~16 near the SuperWASP distance limits) 30% and 75% will be lost 
due to blending for b>20° and b<20° respectively. For the early M dwarfs (V~18), ~50% 
and 85% will be blended for b>20° and b<20° respectively (see Section 3.4). We therefore 
chose to consider only regions of sky with b>20°, and combined this sky coverage with our 
distance sensitivities to derived space volumes. We then combine these volumes with the 
nearby luminosity function of Reid Gizis & Hawley (2002) and the vertical disk density law 
from Zheng et al. (2001), to derived the number of stars within the transit detection limits 
as a function of My. We then factored in a 1% BD companion fraction at a separation of 
0.02AU for late K and M dwarfs (My >8), with a negligible binary fraction for earlier stars 
(see Section 2.1), and thus determined the fraction of binaries that should be fully eclipsing 
using the equations from Section 2.2. 

In order to factor in how many of these BDs will be detectable from the secondary 
eclipse, we simulated the BD Teff distribution. This was done by defining an array of 
~8000 BDs (0. 015-0. 075M0) where we assumed a population mass function with q;=0.5, 
and a uniformly distributed random age distribution from 0-lOGyrs. For each BD we then 
used the NextGen and DUSTY (Chabrier et al. 2000) model isochrones to estimate Tg// from 
mass and age, using linear interpolation between isochrone grid points. We chose NextGen 
Te//s higher than 2500K, and DUSTY Te//S where the NextGen Te//S were <2500K. The 
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resulting histogram of Te//S is shown in Figure 7. Spectral type ranges (some shown on 
the figure) were determined using the spectral type - Teff relations from Reid et al. (1999; 
M6-L2), and Knapp et al. (2004: for >L2). For each of our selected spectral type ranges 
(M6±l-T4±l) we added up all the BDs in the appropriate Tg// range, and divided by the 
total number, so as to represent the fraction of BDs in that range. These results are tabulated 
in Figure 7. The implications of this Teff distribution will be discussed in the next section. 
Finally, we removed any binaries in which the primary was too bright compared to the ED 
companion (using the My limits determined in Section 4.1). 



6.2. Simulation results 



The simulation predicts that a total of ~1000 eclipsing BD companions could be discov- 
ered by SuperWASP over the full La Palma sky. Of these, secondary eclipse measurements 
should allow the BD brightness to be measured (to 10% accuracy) for ~150 systems. The 
predicted distribution of primary and BD spectral types, and BD mass and age for these 150 
systems is shown in Figure 8. Note that these predicted results scale with the assumed BD 
binary fraction and 1/a, and can thus be easily scaled up or down for different input values. 
In the top left histogram, the favoured primary spectral types are the earlier types, since 
these stars are brighter, and detectable out to greater distances. In the top right histogram, 
there is a strong preference for a BD spectral class of L4±l. This preference results from 
three factors in our simulation. Firstly, the spectral type-Tg// relation we used (Knapp et 
al. 2004) gives a relatively large Tg// range for the L3-5 spectral classes (1500-1900K; see 
Figure 7). Secondly, L4 companions will be detectable around all the primaries that our 
simulation considers as potential hosts for BD companions (My >8; see Figure 6). Thirdly, 
although we assume an approximate BD upper mass limit of O.O75M0, the Dusty models 
predict a irhemm ~0.07. Our simulation thus has the higher mass (O.O7-O.O75M0) objects 
stabilising their Te//s around 1500-1900K (ie. L4±l) in the 1-lOGyr age range. Uncertain- 
ties in the spectral type-Te// relation and the precise value of raHBMM could thus change 
the spectral type and width of the peak somewhat. However, note that we certainly ex- 
pect an L4 predominance over later companions, since the L4s are detectable around the 
brighter primaries. The bottom left histogram shows this preference for 0. 070-0. 075M0 ob- 
jects. Finally, the bottom right histogram shows a strong preference for young BDs, because 
young BDs are brighter and more easily detectable from their secondary eclipses. The age 
histogram has ~20 binaries younger than IGyr, ~10 with ages l-2Gyrs, and a continuum of 
~2/Gyr for 2-lOGyrs. This relatively flat age continuum is mainly populated by the L4±l 
type companions. 
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But how many of these detectable BDs would be significantly irradiated? It can be seen 
by looking at Figure 2 that no BDs with spectral types <L0 would be significantly irradiated 
by any of the primaries we consider. However, an L4 BD would be slightly (<25%) irradiated 
by the Mv=8-9 primaries. Later BDs would be irradiated by steadily fainter primaries, and 
the small number of T5 BDs would be irradiated by the majority of primaries. We estimate, 
using Figure 2 and the top two histograms in Figure 8, that ~70% of the predicted companion 
BDs will not be irradiated by more than 10%. 



6.3. A subset of young systems 

We also wish to consider the properties of a young subset of the BDs. Such systems are 
preferred by our detection criteria because BDs are brighter when young. As we discussed 
in Section 5, many young systems should be members of kinematic moving groups, and as 
such have ages that can be well constrained. We therefore re-ran our simulation using a 
Tgff distribution derived by only selecting BDs with ages <lGyr. This simulation predicted 
~20 binary systems in which BDs could be well measured from the secondary eclipse. The 
histograms describing these binaries are shown in Figure 9. The primary spectral type 
histogram is similar in shape to that of the main simulation. The BD spectral type histogram 
however, is somewhat different to the main simulation. The L4 peak is still there, but has 
been significantly reduced in number. This is because although the L3-5 Teff range is 
still relatively large, the majority of the 0. 070-0. OTSM© companions (ie. in the 1-lOGyr 
age range) have been excluded in this age limited simulation. However, there is now an 
additional peak for M5-9 BDs. This peak results from the fact that all M5-9 dwarfs must 
be young in order to be BDs, and will thus not be removed by our age constraint. The mass 
histogram still shows a general preference for higher mass BDs, but the removal of the large 
number of older 0. 070-0. OTSM© objects has decreased the highest mass point significantly. 
The constrained age histogram is shown with a finer scale, and the preference for younger 
systems clearly remains. As before, we have used Figure 2 and the top two histograms (this 
time from Figure 9) to estimate irradiation effects, and estimate that ~80% of the predicted 
young BD companions will not be irradiated by more than 10%. The fraction of the ~16 
binaries (containing un-irradiated BD companions) that will be clearly associated with a 
moving group is unclear (see Section 5.4). However, it seems likely that we might expect 
~10 systems capable of providing BD mass, age, radius, metallicity and luminosity data, for 
our assumed 1% BD companion fraction. 
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7. Identification and followup 

7.1. Detection rate 

On-line optical survey data suggests that we would expect ~300, 000-1 million sources 
to V~18 in a 320 sq. deg area of sky with \b\ > 20° (ie. a typical Super WASP-I pointing). 
Typically ~70% of these will have V=16-18, and ~20% V=14-16. About 20% of sources 
will be late K and early M dwarfs (~100,000 per 320 sq. degs), with the majority being 
earlier G and K dwarfs. Of these, ~20,000 will be monitored with a sufficient accuracy to 
detect BD transits (ie. mostly in the V= 14-16 range). Typically ~30% (~6000) of these 
are expected to be slightly blended (at the >1% flux level) with neighbouring sources. We 
do not expect random fluctuations to produce more than ~1 false detection of a 5-a transit 
event in every ~2 million sources (ie. in the full sky late K and M dwarf population). 
False positives will more likely result from unresolved eclipsing stellar systems in which the 
eclipsing binary component contributes only a small fraction of the total light. Alonso et 
al. (2004) report the discovery of 1 transiting planet out of 16 candidates (where the false 
positives were dominated by such sources). 

We might expect ~5 genuine fully eclipsing BD companions to be discovered in each 
~320 sq. deg area of sky (or for instance, 30 per year from Super WASP-I) . Of these, we 
would expect ~15% of binaries (~5 per year) to have measurable (10% accuracy) secondary 
eclipses, and ~l/3 of these to have ages <lGyr (ie. 1-2 per year). These BD detection rates 
may be scaled for alternative binary criteria - they will scale with BD companion fraction 
and 1/a. Note also that the actual transit candidate numbers will be higher, since we have 
only considered full eclipses, and also expect additional transits by faint low-mass stellar 
companions, and possibly planetary companions to pass our selection criteria. However, the 
level at which these will occur is beyond the scope of this paper. 

7.2. Followup procedures 

The task of identifying eclipsing BD companions will begin with a search of the transit 
survey database. Automated search algorithms should be used to identify transit signatures 
in the light curves, allowing for the possibility that there could also be some variation in 
brightness due to rotationally modulated spot coverage. Large pixels may result in the 
blending of some sources, so it will be important to confirm the transit source at better 
spatial resolution. Optical (eg. from the STSci Digitized Sky Survey, or the Super Cosmos 
Sky Survey) and near infrared (from 2MASS) sky survey images will show if a transit source 
is a blend of several sources or a single star, and also provide optical/NIR colours and basic 
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astrometry. These data should also differentiate between dwarfs and giants, confirming 
or otherwise the faint nature of a transiting companion. It may well also be desirable at 
this stage to measure additional optical light curve data (using a larger telescope), to both 
confirm the transit source, and more accurately measure the primary eclipse profile. With 
the addition of a measured spectral type (and metallicity) , the character of the transit will 
have been well established. 

Having confirmed a low luminosity eclipsing companion of a cool dwarf, the next step 
will be to attempt to measure the secondary eclipse profile and obtain high resolution radial 
velocity data (see Sections 4.1 and 4.2). A parallax distance to the primary will also be 
required to confirm its brightness (from which the BD brightness is inferred). Together these 
data could provide BD mass, radius, metallicity and brightness. The final stage will be to 
try and determine the age of the system, by measuring accurate astrometry and considering 
all relevant moving group membership criteria (see Section 5.3). 



7.3. Additional issues for young sources 

Many young stars may show periodic photometric variability resulting from rotationally 
modulated spot coverage. This generally occurs at the 5-10% level (Krishnamurthi etal 1998), 
has periods of ~l-5 days (Terndrup et al 2000), and maintains the same modulation shape 
over numerous rotation periods (Young et al. 1990). The binary systems we consider should 
produce primary eclipses at least 4% deep, lasting for ~1.5hrs, with a period of ~1.5d. With 
such distinct forms of variation, automated search algorithms should be able to identify 
transit sources even if they show some spot modulation. 

Accurately measuring the much fainter secondary eclipses from foUowup NIR light curves 
of variable sources could prove more of a challenge. If the primary rotation and BD orbital 
periods are different, one could simply determine the form of the rotational modulation by 
phasing the data appropriately, and averaging over several rotation periods. The rotational 
modulation could then be subtracted off the light curves, leaving only the secondary eclipse 
dips. However, if the rotation and orbital periods are the same (close binaries can become 
phase locked by tidal forces; see below), then the eclipse times will always correspond to the 
same rotational phase. 

For a O.O5M0 BD at 0.02AU separation from a primary with mass from O.G-O.SSMq 
(late K - early M), the spin-orbit synchronisation time-scale varies from ~0.04-2.5Gyrs 
respectively (depending primarily on the primary radius and companion orbital separation; 
Udry et al. 2002). Such close BD companions around late K dwarfs are thus likely to be 
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in synchronised orbits, but those around early M dwarfs will probably not have had time 
to synchronise. To deal with synchronised sources, one could measure simultaneous R-, I- 
and K-band light curves. The optical bands would contain a negligible fraction of BD flux, 
and any variability would thus be purely due to rotational modulation. The colour-colour 
relationship between R-K and R-1 modulation away from eclipse could be well determined 
from the light curve data. This relationship will depend on how the flux in each band 
changes as spots rotate onto, around and off the observable face of the star. This will of 
course depend on the temperature, size and surface distribution of spots, but should be a 
tight relationship for any particular spot configuration. This colour-colour relationship could 
thus be used to convert R-I measurements (taken during eclipse) into R-K measurements (as 
if the BD was not eclipsing). The R-band light curve could then be converted into a "pseudo- 
K-band light curve" (using this R-K colour), which would show only the K-band rotational 
modulation during eclipse. This could then be subtracted off the original light curve as 
before. Simultaneous optical/NlR observations would be best done with a robotic telescope 
with infrared/optical capabilities. The 2-m robotic Liverpool Telescope for example could 
measure an appropriate set of R- I- and K-band magnitudes in ~25 minutes. This would 
then require light curve coverage of more transits to obtain the required accuracy. 

The radial velocity data could also be affected by rotationally modulated spot coverage. 
As photospheric dark spots and brighter regions rotate over the surface of the star, they can 
produce radial velocity modulation. These rotational radial velocity variations will have an 
amplitude up to vsini, a period equal to the rotation period (Prot), and a oscillatory shape 
that depends on the distribution of spots across the stellar surface. For illustration, early 
M stars with ages <100Myrs fall into three categories of rotation rate. The majority (~60- 
70%) have slower vsini <20kms^^ rotation velocities (with rotation periods of ~l-5days). 
~10-20% arc faster rotators (vsini =30-60 kms^^), and only 5-20% are ultra-fast rotators, 
with vsinz >100kms^^ (sec Mochnacki ct al. 2002, and Tcrndrup ct al, 2000, for young field 
stars and Pleiades stars respectively). The radial velocity oscillation we expect close BD 
companions to induce is ~ ±20kms~^, with a period of ~1.5days. Clearly in the majority 
of cases, the rotational velocity modulation will be smaller than the orbital modulation 
induced by the BD. In general, our task will simply be to deconvolve the two radial velocity 
components. This should generally be easy. For synchronised orbits, one could measure 
additional radial velocity data at a later date, when the spot modulation has changed. 
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8. Conclusions 

In order to assess the prospects of discovering eclipsing BD companions capable of 
yielding mass, age, radius, metallicity and luminosity data, we have assessed what is currently 
known about such systems. Although close BD companions are fairly rare around late K and 
early-mid M dwarfs, there could be a ~1% companion fraction at very close (~0.01-0.02AU) 
separation. Our knowledge of this companion fraction is poor, but if correct, it suggests that 
although the BD desert extends into this regime, it is not as dry as for solar type stars. 

Using a simulation, and the SuperWASP project as an example, we have predicted the 
numbers and properties of close BD companions that could be identified by ultra-wide-field 
transit surveys in the next few years. We assumed a 1% close BD companion fraction at 
0.02AU, a uniform age range from 0-lOGyrs, and a BD mass function of a— 0.5. We estimate 
SuperWASP to be sensitivite to transits around ~1 million late K/early-mid M dwarfs. 
Based on our current knowledge of BD properties our simulation predicts that ~400 BDs 
in eclipsing binaries could be found. It predicts that ~60 of these binaries should have 
secondary eclipses that are measurable from followup NIR light curves, and should thus 
provide direct measurement of the BD brightness. The most common detectable BDs are 
expected to be higher mass BDs orbiting late K-M3 primaries. The majority (~70%) of the 
detectable BDs should not be irradiated by >10%. 

The simulation also reveals a preference for youth amongst binaries containing detectable 
BDs. This results from the fact that younger BDs are brighter relative to their primary. 
Our simulation predicts that ~20 binaries containing detectable BDs will have ages <lGyr. 
We suggest that it is likely that a significant fraction of such young systems could be young 
moving group members - disk populations with characteristic space motions that are believed 
to share a common origin, and thus have an age that can be well constrained (see Section 
5). The typical spectral type of the BDs in this young sub-sample ranges from M5-L5, and 
may be bi-modal, peaking at M5-9 and L3-5. Higher mass BDs are still preferred, but the 
preference is not as strong for the young subset, and the population is more evenly spread 
over the mass range 0. 040-0. OTSM©. We expect ~80% of these young companion BDs not 
to be irradiated by >10%. 

Although some of our input assumptions are rather uncertain, we show that transit 
surveys such as SuperWASP could find significant numbers of eclipsing BD companions, and 
by carrying out the followup measurements we describe, it would be possible to accurately 
determine mass, age, radius, metallicity, and luminosity data for ~10-20 young (<lGyr age) 
BDs. The majority (~80%) of these should not be irradiated above the 10% level, and could 
thus provide an empirical testbed to compare to theoretical models of free-floating BDs. 
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Fig. 1. — The fraction of companion BDs (O.IR©) that will fully eclipse their primary star 
at a separation of 0.02AU, against primary type. 
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Fig. 2. — Primary My against the companion BD spectral type. The sohd hne represents 
binary systems in which the BDs are irradiated by 10% at a separation of 0.02AU. The regions 
above and below the line represent lower and higher levels of BD irradiation respectively. 
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Fig. 5. — SuperWASP magnitude limits for 5-0" BD (O.IR©) transit detections against pri- 
mary type. Limits are shown for dark (D), grey (G) and bright (B) skies. Sohd hnes 
show magnitudes in the SuperWASP system. Dashed hnes signify the SuperWASP hmits 
converted into V magnitudes. A point from Horne (2003) is shown for comparison. 
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Fig. 6. — Detection and follow-up limits for transiting BD companions. The solid lines show 
the distance limits for 5-a BD transit detections (for dark, grey and bright skies) against 
primary type. The dashed lines indicate distances out to which the mass of a O.OSM© BD at 
0.02AU separation could be accurately measured with V- and I-band echelle spectroscopy on 
an 8-m telescope. The My lower limits shown indicate the maximum brightness a primary 
star may have if we are to be able to measure the brightness of a BD companion (with spectral 
type of T4-M6) with 10% accuracy, from a near infrared hght curve of the secondary echpse. 
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Fig. 8. — Predicted numbers of fully eclipsing brown dwarf companions with measurable sec- 
ondary eclipses, that could be found by SuperWASP. The four sub-panels show the primary 
and BD spectral type distributions, and the BD mass and age distributions. 
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Fig. 9. — The same as figure 8, but for systems witfi ages <lGyr. 
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Table 1: Brown dwarf companion fractions as a function of primary type. Companion frac- 
tions have been estimated from the following; (1) Halbwachs et al. (2000); (2) Duquennoy 
& Mayor (1991); (3) Vogt et al. (2002); (4) Blundell et al. (2004); (5) Els et al. (2001); 
(6) Lui et al. (2002); (7) Gizis et al. (2000); (8) Wilson et al. (2001); (9) Reid & Mahoney 
(2000); (10) Marchal (2003); (11) Nidever et al. (2003); (12) Oppenheimer et al. (2001); 
(13) Reiners (2004); (14) Pinfield et al. (2003); (15) Burgasser et al. (2003); (16) Close et 
al. (2001); (17) Reid et al. (2001). 
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Table 2: Summary of the kinematic moving groups that make up the local young disk pop- 
ulation. 



Moving group 


U, V, w 


Age 
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Table 2: Continued. 
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NGC 1901 cluster 
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